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For a planar FEL configuration we study stimulated coherent spontaneous emission driven by a 
gradient of the bunch current in the presence of different levels of noise in bunches. To perform a 
vast amount of simulations required for obtaining statistically valid results, we developed a memory 
and time efficient one-dimensional simulation code based on the integral solution to a Klein-Gordon 
equation describing the field evolution. The longitudinal granularity of the electron bunch density 
originating from shot noise is maintained throughout the analysis. Three-dimensional effects like 
transverse emittance and diffraction are taken into account in simulations via an effective FEL 
parameter calculated from Xie's fitting formula. Calculations are performed for an FEL model with 
the SwissFEL injector bunch parameters. It turns out that a reduction of noise by several orders 
of magnitude below the level of shot noise is required to mitigate the noise effect. We propose a 
novel scheme that allows for formation of electron bunches with a reduced level of noise and a high 
gradient of the current at the bunch tail to enhance coherent spontaneous emission. The presented 
scheme uses effects of noise reduction and controlled microbunching instability and consists of a 
laser heater, a shot noise suppression section as well as a bunch compressor. The noise factor and 
microbunching gain with and without laser heater are estimated. We found that shot noise reduction 
by three orders of magnitude can be achieved for a finite transverse size electron bunch. 

PACS numbers: 41.60. Cr, 29.27.-a, 41.85.-p, 41.60.Dk 
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I. INTRODUCTION 

Operating free electron lasers in the XUV and X- 
ray regions are usually based on the SASE (Self Ampli- 
fied Spontaneous Emission) process starting from spon- 
taneous emission of undulator radiation initiated by shot 
noise [1 j. The latter is a stochastic process and the SASE 
FEL produces a series of random superradiant spikes 
with a large variation of intensities [2]. A method of 
direct FEL seeding by a coherent quantum laser operat- 
ing in higher harmonic generation regime (HHG) as well 
as methods of high-gain and echo-enabled harmonic gen- 
eration (HGHG and EEHG) in FELs have been proposed 
to overcome this deficiency [3}{5]. Apart from these ex- 
ternal seeding, an intrinsic seeding originating from co- 
herent spontaneous emission (CSE) driven by the cur- 
rent gradient of electron bunches can be employed [6H2] • 
However, shot noise competes with the coherent seeding 
and reduces the radiation pulses coherence and pulse to 
pulse reproducibility. To aid the FEL process recent pa- 
pers |T0l [11] have proposed schemes to decrease the noise 
level below the shot noise level resulting in the so-called 
'quiet' bunches. Specifically, a space-charge dominated 
interaction region followed by a dispersion section sup- 
presses shot noise across a wide range of frequencies (the 
experiment [12] demonstrated that after 1/4 plasma os- 
cillation period the initial periodic THz modulation is 
completely washed out and one obtains a uniform tem- 
poral profile of the bunch density as predicted in [10 J. 
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Recently, the shot noise suppression technique in the op- 
tical wave range was verified experimentally for the first 
time [T3] . 

In view of the shot noise reduction schemes, the in- 
trinsic seeding driven by the current gradient of 'quiet' 
electron bunches becomes an attractive way of obtain- 
ing pulse to pulse stable stimulated CSE radiation with- 
out external seeding lasers. The temporal structure of 
the radiation pulses can be predetermined by prepar- 
ing electron bunches with a special current distribution 
so that amplified radiation pulses will exhibit a repro- 
ducibility determined by that of the electron bunch cur- 
rent [14]. For example, in the experiment [15] very stable 
from pulse to pulse radiation signals very observed from 
rectangular bunches with steep rise and fall of the parti- 
cle distributions at the beginning and at the end of the 
bunch. Thus, our goal is to investigate the stimulated 
CSE in the presence of different levels of noise and de- 
termine the level of reduction of shot noise required to 
obtain radiation pulses with a well predetermined tem- 
poral structure. 

In order to initiate stimulated CSE electron bunch cur- 
rent gradients have to be sufficiently large, this implies 
that the particle distribution at the tail of bunches has 
to have a steep fall. To create bunches with such a distri- 
bution, we propose a novel bunch formation scheme that 
uses the effects of noise suppression and controlled mi- 
crobunching instability p~6j [17] . The latter is employed 
in a way similar to that used in the longitudinal space 
charge amplifier (LSCA) [18 . The principal difference 
from an LSCA will be that only the tail of the bunch 
is allowed to be sensitive to the microbunching instabil- 
ity whereas the main core of the bunch has to be stable 
against the instability. This can be realized by using a 
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laser heater [19] [20] with a partial overlap between elec- 
tron and laser pulses such that the bunch tail remains 
unheated. Then, after a magnetic compressor the bunch 
will have the microbunched tail with a scale of the cur- 
rent gradient variation of the bunch core of the order of 
the variation of the laser used for heating. Further, a 
short noise suppression scheme has to be applied to pro- 
duce 'quiet' bunches with high current gradients. Note 
that magnetic bunch compression used in [15 resulted 
in steep variations of the bunch particle density and al- 
lowed coherent spontaneous emission to dominate over 
incoherent one. However, a detailed study of the forma- 
tion of electron bunches with special density distribution 
is beyond the scope of the present study and we will limit 
ourselves with order-of-magnitude estimates of the bunch 
properties as the bunch passes through the laser heater, 
shot noise suppression section and bunch compressor. 

CSE from 'quiet' bunches is of interest not only for 
VUV or X-ray FELs but also THz FELs may benefit from 
this approach. Moreover, a THz FEL with moderate rela- 
tivists bunches can be an excellent stand for testing shot 
noise suppression techniques and formation of bunches 
with high current gradients via controlled microbunch- 
ing instability. So, we first develop a simple absolute 
stable method for solving the non-averaged equations of 
free-space and waveguide FELs and then apply our sim- 
ulation technique to study lasing in an FEL seeded by 
'quiet' bunches with high current gradients. We will per- 
form simulations for a set of the electron bunch parame- 
ters obtained at the SwissFEL injector. 



II. INTEGRAL SOLUTION TO THE ID 
KLEIN-GORDON EQUATION 

Using a complete set of transverse modes the problem 
of excitation of electromagnetic radiation by relativistic 
electron bunches moving in an undulator in the presence 
of a waveguide or medium can be reduced to a system of 
equations in the form of a ID Klein-Gordon equation [21]. 
Therefore, we will build an integral solution to the fol- 
lowing equation 
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where z and t are the longitudinal coordinate and time, 
respectively. In ID free-space approximation (k± = 0) 
Eq. is the wave equation and A± is nothing but the 
transverse component of the vector-potential, I± is the 
current density averaged over the bunch cross-section. 
In general case, k± and A± are the transverse wavenum- 
ber and mode amplitude resulting from an expansion of 
the vector-potential into a convenient set of transverse 
modes (waveguide modes, Hermite-Gaussian modes, hy- 
brid optical- waveguide modes); I± is the effective current 
density dictated by a convolution of the physical current 
density with transverse modes. 



We label the discrete pointlike electrons by the index 
q = 1, ...,Q e , where Q e is the total number of elec- 
trons in the bunch. Note that Q e is different for dif- 
ferent bunches because of shot noise. Choosing the axial 
coordinate z to be the independent variable, we denote 
the arrival time of the qXh electron at z by t q (z), and the 
transverse and longitudinal velocities by v±\ q (z) and v z \ q , 
respectively. Then, the effective current density reads 



z \q 



S[t-t q (z)\, 
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where S(z) is the Dirac delta- function. In ID free-space 
approximation the sources of radiation can be treated as 
charged disks of infinite transverse extent, with charge 
per unit area g (g < 0). 

In what follows, we omit the vector notations because 
in linearly or helically polarized undulators the trans- 
verse electromagnetic field can be described by a scalar 
equation of the Klein-Gordon type [2T] , 

In an FEL incoherent undulator radiation from rela- 
tivistic electrons is about 7 4 times larger in the forward 
direction (with respect to the bunch propagation) than 
in the backward direction, where 7 is the average bunch 
energy in terms the rest mass. In the case of stimulated 
emission the forward radiation is also strongly dominant 
and we will take into account only it. Then, a solution to 
Eq. ([I]) describing the field co-propagating with bunches 
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A ± (z,t) = ^ J dt' J dz'G(z- z',t-t')I ± (z',t'), 
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where the Green function satisfies the equation 

d^_^& 2 _ 
k8z 2 c 2 dt 2 
and reads 22 



G 



k 2 ± )G = -4tt5(z - z')5(t - t') 



G(z-z',t-t') = 2ircU[c(t - t') - (z - z')] x 

Jo[Wc 2 (i-t') 2 -(*-*') 2 ]. (4) 

Here T e is the bunch duration (electrons enter the in- 
teraction region during the time interval from — T e to 
0), Jo(t) is the Bessel function of the first kind, U[t] is 
the unit-step function. Note that only radiation emitted 
at time t — (z — z^/v^ by electrons in position z' con- 
tributes to the radiation observed at position z and time 
t because of causality. Here, v p h is the phase velocity of 
the radiation field. 

Using Eqs. Q-Q after some algebra we get the electric 
field of the form 
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E / dz 



c~df 

, v Mq (z') dG{z-z',t-t q {z')] 



v z \ q {z') dt 

U[t-t q (z / )]U[t q (z , )^T e }. (5) 
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From Eq. ([5| in ID free-space approximation for given 
electron trajectories in a planar undulator we recover the 
well-known result [8 



E ± (z,t) = -Y,i 



JC 2 /2 



[uj q {t-t e q - z/c)]x 



U[c(t-t e q )-z}U[z-v zlq (t-t e q )}, (6) 

where uj q — 2^k u v z \ q / '(1 + /C 2 /2); /C and k u are the un- 
dulator parameter and wavenumber, respectively; t e , v z \ q 
and 7 g are the entrance time, average longitudinal veloc- 
ity and energy in units of the rest mass of gth electron, 
respectively. We see that the electric field is a superposi- 
tion of plane waves with different frequencies that depend 
on the electron energy. Equation ([6| provides a simple 
insight into formation of coherent spontaneous emission 
in an FEL and is used later for verification of the numeri- 
cal simulation scheme. In particular, for a monoenergetic 
electron pulse with a rectangular current profile Eq. ([6| 
shows that a sequence of radiation pulses, which are one 
radiation cycle in duration, develops at each end of the 
electron pulse, see Fig. 1 and discussion in [9]. 

A plane wave in (|6| produced by qth. electron at longi- 
tudinal position z is nontrivial only if t G \t e q + z/c, t e + 
z/v q \ z \. This result comes directly from the Green func- 
tion Q (recall that Eq. ([6| is derived for k± = 0) and 
it turns out that two unit-step functions in ([5| lead to 
redundant information. Keeping in mind the last obser- 
vation one can check that these two unit-step functions 
m Eq. ([5]) may be omitted without the loss of generality. 



III. NON-AVERAGED FEL MODEL: 
NUMERICAL ALGORITHM 

Let us introduce the traditional set of dimensionless 
variables 23 
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Gain length i g is related to dimensionless Pierce param- 
eter p as i g = l/(2k u p), and p reads 
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where <j r is the rms transverse size of the electron beam. 
Here j r is the resonant energy in units of the rest mass 



JJ ~ Jo (lfuo) ~ Jl (lTuv) 

is the traditional JJ-factor, v\\ is the average longitudi- 
nal velocity, e and m e are the electron charge and mass, 
respectively. 
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Then, using the normalization (|7|), at the Compton 
limit /)< 1 we get the dimensionless self-consistent sys- 
tem of equations in the form 
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where fi\\ is the number of electrons per unit r enter- 
ing the interaction region at ( = and Xq 1S the charge 
weighting parameter Xq — -Tj_(C = 0, T)//_i_ >p k, where 
I±(C = 0? r ) * s the effective electron current with maxi- 
mum value /j_, p k- Note that deriving Eq. ([9| we ignored 
undulator harmonics since CSE at the third harmonic is 
about two orders of magnitude smaller than at the fun- 
damental frequency (higher harmonics are even smaller) 
since the change in the bunch current over the harmonic 
period is smaller and hence the smaller the source term 
driving CSE [14]. 

To avoid bulky expressions we will demonstrate im- 
plementation of the numerical algorithm for ID free- 
space case (k± = 0), then dG[r — r q ((')X — ( f ]/dr = 
S[r — T q ((') — C + C']- Generalization for k± =^ is 
straightforward. It is also advantageous to express in 
an explicit way that the electric field is co-propagating 
with electrons, to this end we define the complex am- 
plitude by the relation F = Re{Fe zu;r ^/ c_ ^} , where 
uj r = 2cj n 7^/(l + K?)/2 is the resonant frequency that 
corresponds to the case of ideal synchronism of bunches 
with the radiation field at the undulator entrance. Then, 
we can re-write Eq. (JoT) as 



dr. 



q = 2pp q , ^ = -2Re{F(C, r q )e^/^ (l - e */') 



(10) 



9=1; 



Note that our system (|10| is equivalent to the correspond- 
ing one of Krinsky [8 and differs from that of McNeil et 
al. [9] by term (l — e^/ p ), which is omitted in [9] and 
relevant in the case of short bunches [8]. 

For the numerical simulation we introduce a uniform 
mesh over £ with fixed step A£ such that the kth node is 
located at = k A(. Then, we replace the integral from 
to £ by a sum of integrals over small intervals A( 

K (fe+l)AC 

/ / (...)dC, (ii) 



k AC 



and expand r q (£') into a Taylor series around k A( up to 
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the linear term 



Adams-Bashforth-Moulton method 



r q (C) *r q (k AC) + ^ 



(C'-fcAC) = 

C'=fcAC 

f)+2 W ( fe )(C'-fcAC) ) (12) 



where r q = T q (kA£) and /4 = /i q (kA£). Then, the 
field amplitude at Kth node may be written as 
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o [r + C' - K AC - rf > - 2^) (C' - fc AC)] dC'. (13) 



Because of the Dirac Delta function, the integral over (J 
is nontrivial only if 



k A( <C(r) < (fe + l) AC, 



where 
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2/9/Xg^ - 1 
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Note that C' depends on time r and we have to scan over 
r to calculate the integral as a function time. Hence, we 
introduce a mesh w.r.t r with step Ar and label nodes by 
subscript j such that tj = j Ar. When inequalities (14) 



are fulfilled the field amplitude in position Cfc at time re- 
reads 
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Note that F(Ck+i, Tj) and F((x^j) are related by 
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As we mentioned one needs to scan the condition ( [l4| ) 
over time to get the field amplitude as a function of time. 
Then, the question is what the time interval of scan- 
ning should be. Electrons enter the interaction region 
from — T e to 0, then assuming the constant velocity one 
finds that the radiation field is nonzero at position z if 
^puise(^) £ [z/c — T e , z/vA. In practice, this condition is 
exact since the tail of the bunch moves with constant av- 
erage velocity v\\ (the tail experiences no radiation field 
because of the slippage) but the head starts to radiate 
field that propagates with the velocity of light. It is 
also clear that electron transit time t q is within inter- 
val tpuise(^) an d we need to calculate field amplitude F 
only for interval t pu \ se (z) at given z. The corresponding 
interval for dimensionless time is dictated by r m i n = 
and r max = C + T e . Thus, the developed technique has 
the advantage over finite-difference techniques [9, 25[ [26] 
that the radiation field is calculated only over the time 
interval where it is nonzero, whereas the use of the finite- 



difference method for solving the excitation equation ( 10) 



would require to calculate the field in all points of a mesh 
over interval r G [0, Cmax + ^e] 7 thus increasing the calcu- 
lation time and required computer memory. Moreover, if 
a full time-space structure of the field is not needed, one 
can calculate the field only on a time mesh surrounding 
the bunch, thus even more decreasing the computational 
time. 



IV. STIMULATED COHERENT 
SPONTANEOUS EMISSION 
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The last formula plays a critical role for fast numerical 
simulations without limitations on the computer memory 
since it allows one to keep in the computer memory only 
2Q e positions of electrons and the field generated at the 
previous step over (. 

The equations of motion are integrated using a com- 
bination of the predictor-corrector scheme and two-point 



We will employ a model of the bunch with a rectan- 
gular shape of the current profile and a uniform random 
distribution of electrons. The rectangular current distri- 
bution may be used in simulations if the rise and fall of 
the bunch density occur at a scale much smaller than 
the FEL resonant wavelength that can be a case for the 
SwissFEL injector if the controlled microbunching is ap- 
plied. To model shot noise we use the technique proposed 
in [27] that put the real electron ensemble into a corre- 
spondence to a smaller ensemble of macroparticles with 
the same statistical proprieties. In our simulations we 
employ a ID non- aver aged FEL model, however, the 3D 
effects are taken into account by using an effective value 
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of the FEL parameter p calculated with Xie's fitting for- 
mula [28 . This formula accounts for the effects of emit- 
tance ('matched' bunch focusing is assumed) and energy 
spread as well as for diffraction. The later turns out to 
be the most severe factor of degradation of the FEL pro- 
cess such that the effective 3D FEL parameter, p e /, is 
almost two times smaller than the ID FEL parameter, p. 
Since the energy spread is accounted for by p e f only sev- 
eral thousands of macroparticles are used in simulations. 
This fact along with our memory efficient mathematical 
algorithm results in fast simulations. Main SwissFEL in- 
jector and proposed FEL parameters are listed in Tables I 
and II [29]. 

Table I. SwissFEL injector parameters 
after bunch compression. 



x 10 



Parameter 


Symbol 


Value 


Electron energy 


j r m e c 2 


250 MeV 


Bunch charge 




200 pC 


Bunch peak current 


h 


350 A 


Transverse rms size 




55 pm 


Energy spread 


a^rrieC 2 


64 keV 


Normalized emittance 


£n 


0.36 mm- mrad 


Bunch length 


T h 


190 fsec 



Table II. Main parameters of the FEL. 



Parameter 


Symbol 


Value 


Undulator period 




4 cm 


Undulator parameter 


1C 


3.2 


Number of undulator periods 


N u 


200 


Normalized interaction length 


Cm ax 


23.83 


FEL wavelength 


A r 


0.511 pm 


FEL parameter 


9 


0.0095 


Cooperation length 


L c 


4.29 pm 


Gain length 


* B 


33.57 cm 


Normalized bunch length 


n 


13.3 


Effective number of electrons 


n \\ 


1.25 -10 9 


Effective energy spread 


p 1 


0.016 


Effective emittance 


p £ 


0.0025 


Diffraction parameter 




2.6 


Effective FEL parameter 


Pef 


0.0053 



Coherent spontaneous radiation emitted by the bunch 
tail is amplified by electrons as it propagates through the 
bunch resulting in the stimulated CSE, see Figs. [l]and|2| 
In Fig. [l] we presented the radiation power as a function 
of time for an idealized case of the FEL without shot 
noise in the electron bunch. In the nonlinear regime the 
radiation pulse has a clear spike with a duration of the 
order of the cooperation length, Ar ~ 1 (At ~ L c /v\\). 
This main spike leaves the bunch behind approximately 
in the cross-section C ~ T e + 1 (z ~ v\\T e + L c ) and then 
it propagates without distortions in free-space. 

Let us remark that in a THz FEL the cooperation 
length may achieve values smaller than 10A r such that 
the main spike of stimulated SCE will be around ten cy- 
cle in duration or even shorter, therefore quite short (as 
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FIG. 1: Normalized power vs. normalized time at different 
cross-sections of the undulator. Shot noise is ignored in this 
simulation. 
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FIG. 2: Typical dependence of the normalized power on the 
normalized time at different cross-sections of the undulator 
with taking into account shot noise. 



compared to the resonant wavelength) pulses can be ob- 
tained by means of the FEL process itself. 

Shot noise results in substantival distortions of the ra- 
diation pulse shape and reduction of the emitted power, 
see Fig. |2j To mitigate the shot noise effect several tech- 
niques have been proposed. However, to the authors 
knowledge the problem to what extent the bunch noise 
has to be suppressed to obtain well predetermined radi- 
ation pulses in the FEL was not studied and we address 
this question to our simulations of stimulated SCE. 

To characterize the level of noise at wavenumber fc, we 
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will follow the ref. [IT] and define the noise factor as 

r (M = -^-J2 eik[Zq(s) ~ Zp(s)] ' ( 19 ) 

where z q (s) is the longitudinal bunch coordinate of par- 
ticle q at position s in the undulator. One can check 
that if the particle positions are uncorrelated, then shot 
noise results in (r(fc, s)) = 1, where (. . .) stands for the 
statistical averaging. The situation (r(fc, s)) < 1 corre- 
sponds to the case of anticorrelated ('quiet') bunches. If 
the bunch energy spread is small, then magnetic com- 
pression causes the microbunching instability so that 
(r(fc, s)) ~ G, where G is the microbunching gain. This 
situation corresponds to the correlated beams case. In 
what follows, we will consider not only 'quiet' bunches 
but also correlated ones since even in the presence of a 
laser heater the microbunching gain can greatly exceed 
unity at certain wavelength when a multistage compres- 
sion is used. 

To study the effect of noise on the FEL per- 
formance, we calculated the density probability dis- 
tribution of the maximal power of the pulse, see 
Fig. [3j As one can expect, the mean maximal power, 



M(|F max | 2 ), increases and the relative standard devia- 
tion, cr(|F max | 2 )/M(|F max | 2 ), decreases as noise reduces, 
see Fig. [4j The shape of the density probability distribu- 
tion also changes significantly. Noise with the level corre- 
sponding to shot noise almost halves the mean maximal 
power as compared to the case of zero noise level (7.8 
of |^max| 2 against 11.7 of |F m ax| 2 )- In the latter case 
the maximal dimensional power is around 5.5 GW and 
the energy stored in the main radiation spike is around 
100 /iJ. For (r(fe)) = 10 -3 the relative standard devia- 
tion of the maximal power from its value for (T(k)) = 
is less than 5%. It turns out that for (r(fc)) <C 1 
the logarithm of the standard deviation, (J (\F max \ 2 ) , 
demonstrates quite linear dependence on (T(k)) so that 
^d^maxl 2 ) is proportional to (T(k)) 1 ^ 3 . We see that 
0"(|^max| 2 ) diminishes slowly as (T(k)) decreases and the 
reduction of the noise level by several orders of magni- 
tude is required to obtain the relative standard deviation 
of the maximal power, let say, less than 10%. A typi- 
cal fourier transform of the normalized electric field for 
(T(k)) = and (T(k)) = 1(T 3 is presented in Fig. [HJ 
Analysis of numerous fourier transforms of F(r) indi- 
cates that noise almost does not affect a lower frequency 
part of the spectrum but can distort a higher frequency 
part. 
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FIG. 3: Histograms of the probability density distribution 
of the normalized maximal power \F\ 2 for different levels of 
noise. Calculations have been performed with 1600 indepen- 
dent statistical events. The bunch and FEL parameters are 
given in Tables I and II. 



V. SHOT NOISE SUPPRESSION AND 
CONTROLLED MICROBUNCHING 
INSTABILITY 

In this section we propose and analyse a scheme that 
allows for formation of 'quite' electron bunches with a 
high gradient of the current at the bunch tail. Let us con- 
sider the evolution of a 'noisy' electron bunch from the 
injector exit to the undulator entrance allowing the bunch 
to pass through the laser heater, shot noise suppression 
section and bunch compressor, see Fig. [6] While the pre- 
vious section describes the FEL performance in a rigorous 
way within ID approximation, in this section we will give 
order-of-magnitude estimates of the relevant bunch char- 



— 10 



10 



cc 10 




10 10 



10 

Noise Factor <r(k)> 



10 



10 



FIG. 4: The mean and standard deviation of \F\' 
noise factor. Both axes are in a logarithmic scale. 
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FIG. 5: Absolute value of the Fourier transform of F(t) 
versus the frequency normalized to the resonant frequency, 

UJ r . 



acteristics during its formation in the proposed scheme 
since our intention to understand the physics of the pro- 
cess rather than to present technical specifications. More 
precise results for accurate technical specifications would 
require start-to-end simulations from a photocathode to 
the FEL undulator entrance. For our analysis we will 
adopt typical approximations used for calculations of the 
microbunching gain and laser heater effect [20 . 

The temporal profile of the photocathode drive laser 
exhibits random fluctuations and they cause the density 
modulations of electron bunches. Because of this inhomo- 
geneity of the electron density the space-charge oscilla- 
tions are initiated and the conversion of density modula- 
tions into energy modulations as well as the reconversion 



occur each 1/4 plasma oscillation period. Then, at the 
injector exit electron bunches have some density modu- 
lation that can be characterized by a bunching factor 



b(k) 



1 



Q e ec 



I(s) 



-iki 



ds, 



(20) 



which is related to the noise factor defined, by 
T(k) = Q e \b(k)\ 2 . Here I(s) is the bunch current. The 
microbunching gain is defined as a ratio of the bunching 
factors after and before compression. Following assump- 
tions [20 j, in what follows we neglect the bunch energy 
modulation accumulated inside the injector and analyse 
the evolution of the electron density modulation starting 
from the injector end. 

Let us estimate (T(k)) using Eqs. (13), (22)-(24) from 
the ref. [11] (note that the arguments of the Bessel func- 
tions should be o r kj^ instead of the dimensional argu- 
ment k). The noise factor at the wavelength region from 
A r /2 to 3A r /2 is presented in Fig. [7j Recall that A r 
is the FEL resonant wavelength. Because of the finite 
transverse size of electron bunches, (T(k)) depends on 
the wavelength and we chose dispersive strength R 5 q to 
minimize (T(k)) at A = A r . In our calculations, the drift 
section, Ld, is 10 meters. From Fig. [7] one can say that 
the noise reduction in the scheme employing the disper- 
sive section [TT] is quite uniform as compared to the noise 
reduction obtained only by using the drift section of ap- 
propriate length [30] (cf. our Fig. [7] and Fig. 3 in [30 ). 
In the later case, the noise reduction is observed only for 
a discrete set of wavelength. We conclude that the noise 
reduction by three orders of magnitude is achievable. 
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FIG. 6: The schematic illustration of a possible formation of 

Because of impedances (mainly the space-charge 
impedance) the electron bunch density gives rise to the 
energy modulation as bunches propagate though the drift 
section. This energy modulation is converted into an ad- 
ditional density modulation during magnetic compres- 
sion, so-called microbunching gain. To avoid the mi- 
crobunching gain, the laser heater that increases the un- 
corrected energy spread thus providing strong Landau 
damping against the instability can be used. In Fig. [8] 



quiet' bunch with a high current gradient at the bunch tail. 

we present the results of calculations of the microbunch- 
ing gain for the system with and without laser heater. 
The maximal energy modulation, A7i / m e c 2 , produced 
by the laser heater has to be chosen in such a way that 
it results in an effective microbunching suppression on 
the one hand, and it does not affect the FEL perfor- 
mance because of the increased energy spread on the 
other hand. Without heating of electrons, an appreciable 
microbunching gain is observed approximately at a wave- 
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FIG. 7: Noise factor versus wavelength for a bunch of the 
finite transverse size. 



length of 1 /am. Therefore, if we modulate the bunch tail 
with a period of 1 fim and heat with the laser heater only 
the main core of the bunch leaving the tail unheated, as 
it is shown in Fog. |6j then we will get the microbunched 
tail and the steep variation of the electron density at the 
end of the main bunch core. Note that the microbunch- 
ing gain in the vicinity of the resonant wavelength, which 
is 0.511/im, is less than unity so that the reduced level of 
noise is preserved after magnetic compression. 



Conclusion 

We studied stimulated coherent spontaneous emission 
in a planar FEL driven by electron bunches with a rect- 
angular shape of the charge density and different levels of 
density fluctuations. These density fluctuations originat- 
ing from shot noise compete with an intrinsic coherent 
seeding driven by the current gradient and are in the fo- 
cus of our investigation. We also proposed a scheme that 
allows for formation of electron bunches with a reduced 
level of noise and a high gradient of the current at the 
bunch tail to enhance coherent spontaneous emission. 

The need for a fast non-averaged time-dependent simu- 
lation code able to perform a vast amount of calculations 
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FIG. 8: Microbunching gain versus wavelength with and 
without the laser heater. 



required for obtaining statistically valid results stimu- 
lated us to develop such a code. The mathematical al- 
gorithm behind our code is conceptually similar to that 
behind the code 'Fast' [31] and it employees an integral 
solution via the Green function to the ID Klein-Gordon 
equation governing the electric field evolution. The algo- 
rithm presented is absolutely stable and does not impose 
constraints on discretization steps in longitudinal coor- 
dinate, observation time and entrance time of electrons. 
In the paper the field is calculated only on a mesh sur- 
rounding the bunch and then it is interpolated to the 
electron's positions, thus decreasing the computational 
time as compared to finite-difference methods that uses 
a mesh depending on boundary conditions (advantages 
of integral methods over differential ones for injector de- 
sign codes are discussed in details in [32 ). Our simu- 
lation code was verified against analytical result (|6| in 
a low gain regime and against results presented in [9 
for a deep nonlinear regime and excellent agreement is 
found. To account for 3D effects we used in simulations 
the effective FEL parameter calculated from Xie's fit- 
ting formula. Therefore, our simulation code being ID 
nevertheless takes into account transverse emittance and 
diffraction. 

The structure of electromagnetic pulses in the stud- 
ied FEL is mainly predetermined by the stimulated co- 
herent spontaneous emission since the coherent part of 
spontaneous undulator radiation driven by the current 
gradient is dominant over the incoherent one associated 
with shot noise. However, the incoherent emission be- 
ing small leads to essential distortions of radiation pulses 
during amplification and strongly affects the probability 
density distribution of the maximal power, as shown in 
Fig.[3j It turns out that the standard deviation of |F ma x| 2 
diminishes slowly with the reduction of noise because 
^(l^maxl 2 ) is proportional to (T(k)} 1 ^ 3 for (r(fc)) <C 1. 
We found that shot noise should be reduced by at least 
three orders of magnitude to provide the relative stan- 
dard deviation of the maximal output power less than 
5%. The FEL process starts from spontaneous undula- 
tor emission and the ratio of the power of CSE to the 
power of incoherent emission is Q e |^(^)| 2 , where F^uo) 
is the form-factor of the bunch [7 j. Therefore, for FELs 
operating in other wavelength regions we may expect the 
standard deviation of the maximal output power to de- 
pend on the noise level in a similar way as we found, if 
the value Q e | J^(o;)| 2 is fixed. Thus, VUV or X-ray FELs 
based on stimulated CSE would require the reduction of 
shot noise by three or four orders of magnitude to have 
radiation pulses with a well predetermined time struc- 
ture. 

CSE is driven by the gradient of the bunch current 
and bunches with a steep rise of the electron density at 
the tail are preferable to drive the FEL. To create such 
bunches also having a reduced level of noise we proposed 
a scheme that uses effects of noise reduction and con- 
trolled microbunching instability, and consists of a laser 
heater, a shot noise suppression section as well as a bunch 
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compressor. In our scheme the tail of the bunch is allowed 
to be sensitive to the microbunching instability whereas 
the main core of the bunch has to be stable against the 
instability. This can be realized by using a laser heater 
with a partial overlap between electron and laser pulses 
such that the bunch tail remains unheated. Then, the 
bunch passes through the shot noise suppression scheme 
proposed in [TT] that consists of drift and dispersion sec- 
tions, where the level of noise is reduced. We found that 
shot noise reduction by three orders of magnitude can be 
achieved for electron bunches produced by the SwissFEL 
injector. Then, in the bunch compressor the bunch un- 
dergoes the longitudinal compression and the controlled 
microbunching of its tail occurs as well. We calculated 
the microbunching gain, G, with and without laser heater 
and in the latter case G attains its maximal value around 
80 at a wavelength of 1 fim (start-to-end calculations 
from the photocathode shows that G can even achieve a 
value of several thousands [33]). Therefore, if we mod- 
ulate the bunch tail with a period of 1 jiva by means 
of a photocathode laser (or with an undulator and an 



external laser) and heat with the laser heater only the 
main core of the bunch leaving the tail unheated, as it is 
shown in Fog. |6j then we will get the microbunched tail 
and the steep variation of the electron density at the end 
of the main bunch core. Note that the microbunching 
gain in the vicinity of the resonant wavelength, which is 
0.51 1/xm, is less than unity so that the reduced level of 
noise is preserved after magnetic compression. Thus, one 
can form 'quiet' bunches able effectively to drive CSE in 
the FEL. 
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